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In the present paper, we study spin-^ antiferromagnetic (AF) Heisenberg model 
on layered anisotropic triangular lattice and obtain its phase structure. We use 
the Schwinger bosons for representing spin operators and also coherent-state path 
integral for calculating physical quantities. Finite-temperature properties of the 
system are investigated by means of the numerical Monte-Carlo simulations. Detailed 
phase diagram of the system is obtained by calculating internal energy, specific heat, 
spin correlation functions, etc. There are AF Neel, paramagnetic and spiral states. 
Turning on plaquette term (i.e., the Maxwell term on a lattice) of an emergent U(l) 
gauge field that flips a pair of parallel spin-singlet bonds, we found that there appears 
a phase that is regarded as a deconfined spin-liquid state, though "transition" to 
this phase from the paramagnetic phase is not of second order but a crossover. In 
that phase, the emergent gauge boson is a physical gapless excitation coupled with 
spinons. These results support our previous study on AF Heisenberg model on a 
triangular lattice at vanishing temperature. 

I. INTRODUCTION 

Study of quantum spin models has a long history. In particular after the discovery of 
the high-temperature superconductors, exotic quantum spin states have been intensively 
explored. Among them, spin-liquid state with a deconfined spinon has interested many the- 



oretical and experimental researches [l]. Recently experiments of the anisotropic triangular 
antiferromagnet (AF magnet) CS2CUCI4 revealed the existence of the incommensurate spiral 



order at low temperature (T) and also spinon-like excitations at intermediate T\ 



3|. One 
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may think that this material can be a candidate for so-called Z 2 spin liquid (4 ] . Further- 
more very recently, evidences for a spin liquid in EtMe3Sb[Pd(dmit) 2 ]2 at very low T were 



reported 



5] 



In the previous paper jo], we studied frustrated AF Heisenberg model on an anisotropic 
triangular lattice in two dimensions (2D) at T = 0. We used the Schwinger bosons for 
representing s = | spin operators, and derived an effective model for low-energy region 
assuming existence of a short-range spiral order. Then low-energy excitations are spinons 
and an emergent gauge field with local Z 2 gauge symmetry We studied the effective gauge 
model for the quantum AF Heisenberg model by means of Monte Carlo (MC) simulations, 
and obtained phase diagram. There exist the spiral state, paramagnetic (PM) dimer state 
and the spin-liquid state in the phase diagram. These phases can be also labeled by the 
gauge dynamics, i.e., Higgs, confinement, and Coulomb phases, respectively. 

In the present paper, we shall continue the study in the previous work and investigate 
closely related model, i.e., AF Heisenberg model on layered anisotropic triangular lattice. 
We are interested mostly in finite-T properties of the model. By using the Schwinger bosons 
and CP 1 path-integral methods, direct application of the MC simulations becomes possible 
without assuming any kind of (short-range) orders. We shall clarify the phase diagram of 
the model. Results obtained in this paper support the study in the previous paper [6]. 

This paper is organized as follows. In Sec. II, we shall introduce models and give a 
derivation of the effective gauge-theory model. By using the Schwinger bosons, local U(l) 
gauge symmetry naturally appears. Then we discuss possible phases in the frustrated AF 
spin systems. Section III is devoted to numerical studies. We investigated phase structure 
of the model by calculating the internal energy, specific heat, spin correlation functions, etc. 
These numerical calculations show that there exist AF, PM, spiral and spin-liquid phases. 
Detailed study on the critical behaviors between these phases is given. Section IV is devoted 
to summary and discussion. 

II. QUANTUM AF SPIN MODEL, SCHWINGER BOSONS, CP 1 
REPRESENTATION AND GAUGE THEORY 



In the present paper we shall study a spin- ^ anisotropic AF Heisenberg model on a layered 
triangular lattice shown in FigJTJ For simplicity, we first consider a 3D cubic lattice and 
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FIG. 1: Three-dimensional layered lattice on which the spin models (|2.ip are defined. 

then add diagonal links in the upper-right direction (1-2 direction) in 2D layers. Quantum 
Hamiltonian of the spin system is given as 

H = J\ S x ■ S x+fl + J2 S x ■ S x+ i + 2, (2-1) 

X,fl X 

where S x is spin operator at site x and fi(= 1, 2, 3) is a direction index and also denotes unit 
vector in 3D lattice, whereas 1 and 2 are those of the 2D lattice. Therefore, the Ji-term in 
Eq. (j2.ip is the 3D nearest-neighbor (NN) AF interaction, whereas the J2-term is the next- 
nearest-neighbor (NNN) AF coupling in 2D layers. There exists AF Neel order for J x 3> J2 
at low temperature (T), and it is expected that a quantum phase transition takes place as 
J 2 is increased. 

In this paper, we shall investigate finite- T properties of the system ( 12. ip in detail. To 
this end, we employ the Schwinger-boson representation and the coherent-path-integral 
methods[7j. By means of these methods, numerical study of the system can be performed 
straightforwardly. In terms of the Schwinger bosons at site x, a = (a x ^, a^Y (where O l 
denotes transpose of the vector/matrix O), the spin operator S x is expressed as 

S x = ^a) x aa x (2.2) 

where a is the Pauli spin matrices. As the magnitude of the quantum spin is |, the physical 
states of the Schwinger bosons, \Phys), have to satisfy the following constraint at each site 

x, 

al a a xa \Phys) = \Phys). (2.3) 

<r=t4 

We use the coherent-state path integral for the study of the system (12.11) expressed in 
terms of the Schwinger bosons. To this end, we introduce CP 1 variables z x = (z x ^, z x ^f 
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corresponding to a xa , which satisfy the constraint 



er=t,4. 



as required by Eq. (12.3j) . 

Then the partition function Z is given by 

Z = J [Dz] C p exp drA(r) 
A{r) = -^z X(T d T z x<T -H{z,z), (2.5) 

X,(T 

where r is the imaginary time, /3 = 1/(/ebT), [Dz] denotes the path-integral over CP 1 

variables and H(z,z) is obtained from Eq. fl2.ip by using Eq. fl2.2p . H(z,z) in Eq. fl2.5p is 
explicitly given as follows, 



H(z,z) — — \z x z x+fl \ 2 + — 1^^+1+2 



2 



= -y^|^^+ M | 2 -y2Zkx^+i+2| 2 + constant, (2.6) 

where 5 X = (z x ±, —z x fY, which is nothing but the time-reversal spinor of z x , and we have 
used the fact that z x and z x are an orthogonal and complete set of vectors in the CP 1 space. 
If one tries to numerically study the system (12.51) by means of the MC simulations, one 
immediately encounters difficulties in the important sampling procedure because the first 
term in the action A(t) is pure-imaginary. For J\ ^> J2, it is known that by integrating 
out the CP 1 variables z x at all odd (or even) sites of the cubic lattice by assuming a short- 
range AF order, resultant action has a quartic form of z xa (x G even sites) and has a lower 
bound [8]. This calculation, however, cannot be applicable for the case J\ ~ J2 that we are 
interested in. Therefore we shall take another way to avoid the imaginary term in A(r), i.e., 
we consider finite-T properties of the system (12.51) and ignore the imaginary-time dependence 
of variables z x . Study of finite-T properties of the system is not only interesting itself but 
also gives an important insight into low-T properties of the system as it is expected that 
an ordered phase at finite T survives at lower T's[9|. In this approximation, the partition 
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function is given as 

Z = J [Dz\ c ? exp(-S 



D — 7T- \ Z x z x+fi\ H — \ Z x z x+l+2\ 



2 



2 



—Ci |^^a;+At| 2 + | ^^+1+2 1 2 , (2-7) 



X,fJ, 



where ci = ^ and di = ^f-. 

For later discussion, it is useful to rewrite the action So in Eq. (12.7ft as follows. We first 
rename CP 1 variables at odd site of the cubic lattice as 

z x z x , x G odd site, (2.8) 

and introduce a gauge field U xli at link (x, (i) for the AF spin-pair channel. Then 

Si = -c[ }Xz x+fl U XIJ ,z x + c.c.) + d 1 ) j \z x z x+ i +2 \ 2 . (2.9) 

X,fl X 

One-link integral over the gauge field U XfJ , = e xtl can be performed exactly, 

J exp (c^Za.+nU^Zg. + c.c.j = exp (loglo^l^+^l)), (2.10) 

where Jo is the modified Bessel function. From the behavior of the modified Bessel function 
Jo, relation between the parameters C\ and c[ is obtained as follows, 

, U for Cl »l, , ^ 

ci ~ { (2.11) 

[ (d) 1 / 2 for Cl < 1. 
Action 5*i is invariant under the following local gauge transformation, 

Zxa ^ z xa e , z xa y z X(7 e , U x ^ y c +li U x ^e , (2.12) 

where A x is an arbitrary gauge-transformation parameter. In order to investigate the pos- 
sibility of the appearance of a spin liquid with deconfined spinon excitations, study of the 
gauge dynamics and behavior of the gauge field U Xfl is important and necessary. 
We also add plaquette term of the gauge field U Xfl to the action Si in Eq.( 12.9[) as 

z x z x+ i + 2\ 2 — C2 U x ^U x+ ^ jU U x+Utfl U xu +c.c., (2.13) 

X,/l X X,fl>V 
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where the last plaquette term is the counterpart of the Maxwell term of the gauge field 6 X/1 
in the continuum, and corresponds to the ring-exchange terms of spins like 

(c 1 ) 4 C2(S x ■ S x -|- M )(S :r + v ■ S x +n+i/) + ■ • • , 

for the case of small value of o^. Higher-order terms of Ci correspond to nonlocal interactions 
between spins. It should be noticed that the above plaquette term is defined on the 3D cubic 
lattice, and therefore the induced ring-exchange interaction is three-dimensional. The gauge 
field U Xfl is related to the original Schwinger-boson operators as 

TT a x+ ^a xi - a x+fMi a xt , x G odd site 

Ux » ~ \ *t ~t ~t ~t • ( 2 - U > 

[ a x+rt a k - < +/4 <t' x G even Slte ' 

i.e., U Xfl corresponds to creation and destruction operators of spin-singlet bond at sites x 
and x + fi. Therefore the C2-terms in the action S2 (12.131) flip pairs of parallel nearest- 
neighbor spin-singlet bonds, and enhance appearance of the resonating-valence-bond (RVB) 
liquid 0. 



From the previous studies [a, 111], phase structure of the quantum spin models correspond- 
ing to S2 in Eq. fj2.13j) is expected as follows, 

1. For d\ = Ci = 0, a phase transition from a paramagnetic state to the Neel state with 



AF long-range order takes place as c\ is increased 



12j |. In a gauge- fixed formalism, the 



AF Neel state corresponds to the state in which (z x ) 7^ and {U XfJ ) ~ 1. 

2. As the value of di is increased in the AF phase, a spiral state appears at some critical 
value of G?i c (ci)[13j]. In the spiral state, z x is parameterized as 

z x = ^{e iwx v x + e* M v x ), 

where a; is a constant, and the condensation of smoothly varying field v x takes place, 
(v x ) ^ 0. 

3. Furthermore, as the value of C2 is increased, a spin-liquid state with a deconfined 
spinon appears. In the spin-liquid phase, (z x ) = (v x ) = and the gauge dynamics of 
U XfJj is in the Coulomb phase. Gapless gauge boson appears as a low-energy excitation 
coupled to spinons. 
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In the following section, we shall show the results of study on the phase diagram and physical 
properties of the models S 2 in Eq. fl2.13p and So i n Eq.f l2.7p . which support qualitatively 
the above expectation. As mentioned in the introduction, the experiments for frustrated 
quantum magnet Cs 2 CuCl 4 observed the spiral state and deconfmed spin-liquid state {2, 3|. 
Experimental results suggest a crossover in nature of the excitations from spin-1 spin waves 
at low energies to deconfmed spin- 1/2 spinons at medium to high energies. Furthermore, 
EtMe 3 Sb[Pd(dmit) 2 ]2, which is studied intensively these days, is closely related to the present 
model. Therefore, results in this paper are relevant to these materials. 



III. NUMERICAL STUDIES 



A. C2 = case 



In the previous section, we have derived the effective models of the CP 1 -!- U(l) gauge 
variables from the AF Heisenberg model on layered triangular lattice. In this section we 
show results of the numerical study of the models obtained by means of the MC simulations. 
We employed the free boundary condition in the 1 — 2 plane as the system may have an 
incommensurate spiral order with the layered structure. 

We first consider the case with c 2 = 0. We investigated phase structure of the model S 2 
by calculating the internal energy E and the specific heat C for various values of c[ and d\, 

E = ^(S 2 ), C = ±((S 2 -(S 2 }) 2 }, (3.1) 

where L is the system size of the 3D lattice. In the practical calculation, we employed the 
local update by the standard Metropolis algorithm for the total system with size (2 + L + 
2) x (2 + L + 2) x L and performed measurement of physical quantities in the central LxLxL 
subsystem (l^]. 

We have found that E exhibits no anomalous behaviors, whereas C exhibits singular 
behaviors that indicate existence of second-order phase transitions as c[ and d\ are varied. 
Observed phase transition lines in the d\ — c[ plane are shown in Fig21 

We first consider the PM-AF phase transition. In Figj3j we show E and C as a function 
of d x for d\ = 1.0. It is obvious that E exhibits no anomalous behavior whereas C has a 
peak at c[ ~ 3.75 and the peak develops as the system size is increased. This behavior of 
C indicates a second-order phase transition at c[ ~ 3.75. In order to verify existence of the 
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FIG. 2: Obtained phase diagram for model S2 with C2 = 0. There are three phase transition lines, 
which separate AF, PM and spiral phases. All phase transitions are of second order. Locations of 
phase transition lines are determined by calculations of system size L = 16. 
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FIG. 3: Phase transition from PM to AF phases. E exhibits no anomalous behavior whereas C 
has a peak indicating existence of second-order phase transition. 

phase transition and understand physical meaning of each phase, we investigated correlation 
functions of spins that are given as follows, 

G s{r) = ^ ( Ux ' n *+J>)> G 's( r ) = ( n x ■ ^+(1+2),.), (3.2) 
3=1,2 
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FIG. 4: Spin correlation functions for d\ = 1.0. 

where n x = (z x az x ). Numerically obtained results are shown in FigJH At c[ = 3.6, the 
correlation functions have no long-range order (LRO). On the other hand at c' x = 4.0, 4.7, 
they exhibit AF LRO. (Please recall that we have changed variables z x — >■ z x , x G odd site.) 
From this result, we conclude that transition from the PM to AF phases takes place at 
c[ ~ 3.75. 
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FIG. 5: Phase transition from PM to spiral states. Total specific heat C has a peak that develops 
as L is increased. 



We turn to the phase transition from the PM to spiral states as shown in Figj5j For 
d\ = 2.0, calculation of the total specific heat C as a function of c[ is shown in Figj5] We 
also measured the specific heat of each term of the action, which is defined similarly to C 
in Eq.f l3.ip . in order to see the physical meaning of the phase transition. 

1 



C c 



L ? > 
1 



((s c -{s c )Y), 

((S d -(S d )f), 



(3.3) 
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where 



S c — c \ ^ ^ ( z x+^Ux^j, z x + C.C.) 
z x z x+l+2\ ■ 



(3.4) 



See Figj6l From these results, it is obvious that a second-order phase transition from the 
PM to spiral states takes place at c[ ~ 3.9. It is interesting to see that the c^-term of the 
action tends to fluctuate strongly at the phase transition point but the (ii-term does not. 
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FIG. 6: Specific heat of each term in action. Results show that ci-term in action fluctuates strongly 
at phase transition but di-term does not. 

We measured the spin correlations at c[ = 3.8, 4.1 and 4.2. The results are shown in 
FigJ3 It is obvious that at c[ = 3.8 the spin does not have a LRO, whereas at c[ = 4.1, 4.2 
it has a spiral LRO. One should notice, however, that at c[ = 3.8 the spin correlation has a 
short-range spiral order and therefore we call this "phase" a tilted- dimer state, though there 
is no sharp phase boundary between the ordinary PM (at d\ 1) and tilted-dimer state. 
This observation supports our previous study of the AF magnets on anisotropic triangular 
lattice assuming short-range spiral order 6(. We also measured the spin correlation in the 
inter-layer direction in the spiral state, and found that it has an ordinary AF correlation as 
it is expected. It is interesting to see a snapshot of spin configurations in the spiral and AF 
states. See FigJHl 

We calculated the specific heat and spin correlations for various values of the parameters 
c[ and di and have obtained phase transition line that separates the PM and spiral phases. 

Finally let us turn to the spiral- AF phase transition (see Figj9j). We show calculations 
of C and the specific heat of c[ and (ii-terms. See Figsj9] and [TUl It is obvious that the 
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FIG. 7: Spin correlations in PM (tilted-dimer) and spiral states. For c' x = 4.1 and 4.2, there is a 
LRO close to 120°-Neel order. Please remember that direction of odd-site spins has been reversed. 
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FIG. 8: Snapshots of spin configuration in spiral (left) and AF (right) states. In the snapshots, 
direction of spins at all odd sites is inverted. Arrows indicate direction of spins projected into the 
1-2 plane. 

total specific heat C exhibits only very weak anomalous behavior but C c and Cd both show 
sharp peak at c[ ~ 6.7 as the system size is increased. From this result, we conclude that 
the transition from the spiral to AF phases is of second order. 

It is also interesting to see how the spin correlation changes from the spiral to AF phases. 
Results in FigsJTT] show that the spin correlation gradually changes from the spiral order to 
AF order. 

We also numerically studied the original model So in Eq.f l2.7p . and obtained similar results 
to those of model S2. Obtained phase diagram is shown in FigJT2| and spin correlation 
functions in FigJT3l Phase transitions are of second-order and spin correlation functions 
have similar behavior to those of 5*2. 
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FIG. 9: Phase transition from the spiral to AF phases. Total specific heat C for spiral to AF phase 
transition exhibits only weak anomalous behavior at phase transition. 
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FIG. 10: Specific heat of each term for spiral to AF phases transition. 
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FIG. 11: Spin correlation in spiral and AF states. 
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FIG. 12: Phase diagram of model So- There are three phases, AF, PM and spiral phases as in the 
model S2- All phase transition lines are of second order. 
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FIG. 13: Spin correlation functions in model Sq for d\ = 3.0. 

Result for the spiral state obtained in the present subsection obviously means that the 
120°-Neel state is realized in each layer for the case of the isotropic triangular case J = J'. 
This result is in good agreement with the previous study on the AF Heisenberg model on 



isotropic triangular lattice at T = 



151 ] . On the other hand, some of the previous study 



on the anisotropic AF Heisenberg model on 2D triangular lattice at T = suggested the 
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existence of spin-liquid phases 16[]. However the results obtained in this subsection show that 
it does not exist in the present model. In order to address the possibility of the spin-liquid 
phase, we shall study effect of the plaquette term of the emergent gauge field U xft in the 
following subsection. 



B. c 2 > case 



In this subsection, we shall study the effects of the plaquette C2-term in the action S2. 
For the c 2 = case, we found that there exist three phases, i.e., the AF, spiral and PM 
phases. It is expected that the gauge dynamics in both the AF and spiral phases is in the 
Higgs phase as the condensation of the spinon field z x suppresses fluctuations of the gauge 
field U XfJi . Low-energy excitations are gap less spin wave in the both phases. On the other 
hand in the PM phase, it is known that the confinement phase is realized and low-energy 
excitations are bound state of the spinous like a spin-triplet (z x az x ) because of the strong 
fluctuations of U Xfl . 

As explained in Sec. 2, another possible phase in the quantum spin systems is the spin- 
liquid phase. In that phase, there exist no LRO's, whereas low-energy excitation is the 
deconfined spinon z x . This means that in the spin liquid only small fluctuations of the gauge 
field U Xfl are realized and the gauge dynamics is in the Coulomb phase. Knowledge of the 
gauge field theory suggests that such a spin-liquid phase may be realized by turning on the 
plaquette c 2 -term because this term suppresses large fluctuations of the gauge field [ITT] • In 
the previous study on the Z 2 gauge model of the spiral and spin-liquid phases jsj, we found 
that the deconfined spin-liquid phase is realized in the vicinity of the spiral and tilted-dimer 
states. In the present paper, we show the results of study on the U(l) gauge model S2 in 
the c[ — C2 parameter plane with the value of d\ fixed. 

We first show phase diagram obtained by the MC simulations for d\ = 2.0 in FigJT4l 
Reason for choosing this value of d\ is that the spiral and tilted-dimer states appear as the 
value of c[ is varied for c 2 = 0. As shown in Fig{T3J there exists a crossover line emanating 
from the point (c 2 — 1.5, c\ = 0) in the vertical direction. This crossover line separates dense 
and dilute instanton regions, whereas the both regions belong to the confinement phase of 
the U(l) gauge model in 3D. Besides the crossover line, there exist two sharp second-order 
phase transition lines emanating from (c 2 = 0, c[ = 4.2) and (c 2 = 0, c[ = 5.5), respectively. 
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FIG. 14: Phase diagram in C2~d x plane for d\ = 2.0. Spin-liquid phase appears in the vicinity of 
PM (tilted-dimer) and spiral phases. 

As shown in FigIT4"l these are the spiral and AF phase transition lines, respectively. We also 
found another "transition line" emanating from (c2 — 1.05, ~ 2.2), which we identify as 
a crossover to the spin- liquid phase. 

We show the total specific heat C as a function of c' x for C2 = 1.8 and d\ = 2.0 in FigfT5l 
There are two peaks at c[ ~ 1.4 and 1.8, and the second peak at c[ ~ 1.8 develops as the 
system size is increased indicating a second-order phase transition. Calculation of the spin 
correlation given later on shows that it is the phase transition to the spiral state. On the 
other hand, the first peak at c[ ~ 1.4 does not develop as the system size is increased. More 
detailed calculation is shown in FigJT6l There exists small system-size dependence, but we 
think that this size dependence comes from the free-boundary condition that we took for 
the calculation. 

It is useful to see how specific heat of each term behaves. See FigJTTl The specific heat of 
the C2-term C C 2, which is defined similarly to C c and Cj, is a decreasing function of c[ and 
changes its behavior at c[ ~ 1.4. On the other hand, the specific heats of the c^-term and 
eq-term both have peaks at c[ ~ 1.8 and c' x ~ 5.5 and these peaks develop as the system size 
is increased. This result suggests that there is a second-order phase transition at c[ ~ 5.5 
besides at c[ ~ 1.8. 

It is important to see how the spin correlation function behaves and verify properties of 
each phase observed by the measurement of C . Obtained results of the spin correlation for 
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FIG. 15: C as a function of c[ for C2 = 1.8 and d\ = 2.0. 

C2 = 1.8, d\ = 2.0 are shown in FigJTHl At c[ = 1.0 and 1.5, there exists no LRO, whereas 
at c[ = 2.4, 4.0 and 5.0 the LR spiral order appears. Furthermore at c[ = 6.0, the spin 
correlation shows the AF LRO. All the above results verify the phase diagram shown in 
Figm 

In order to investigate the gauge dynamics, it is useful to study instanton (monopole) 
density p, which measures magnitude of topologically nontrivial fluctuations of the gauge 



field U xi j,. p{x) is defined as follows for the gauge field configuration U x ^ = e^*<* [111 Il8 |. 



First we consider the magnetic flux Q x ^ u penetrating plaquette (x, x + /x, x + /i + u, x + v) 

(-4vr < Q Xtlu , < 4vr). (3.5) 
We decompose Q x ,nv into its integer part n x ^ u , which represents the Dirac string (vortex 
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FIG. 17: Specific heat of each term in action for C2 = 1.8 and d\ = 2.0. 
line), and the remaining part Q x ^ u , 

Q x ,fiu = 1nn X)llv + <3> x ^u, (-7T < Q X)f J,u < 7r). (3.6) 
Then instanton density p(x) at the cube around the site x + | + | + |of the dual lattice is 
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FIG. 18: Spin correlation function for various values of c[ 
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FIG. 19: Instanton density as a function of c^. Arrows indicate location of crossover observed by 
calculation of C . 



defined as 



(3.7) 



where e^ u \ is the antisymmetric tensor. 
In FigJT9| we show the calculation of 

X 

As the gauge dynamics is already in the dilute-instanton region of the confinement phase for 
c[ = 0, the value of p is small, but it decreases at c[ ~ 1.4, and its behavior becomes clear 
as the system size is increased. This result indicates that the region between two peaks at 



19 



c[ = 1.4 and 1.8 corresponds to the "deconfined Coulomb phase". In this phase, the gapless 
gauge boson 9 Xfl appears as a low-energy excitation besides the deconfined spinons. 

The global phase diagram in FigJTH is consistent with that of the Z 2 spin-liquid model 
obtained in Ref . 6] . In the Z 2 model, however, the phase transition from spiral phase to PM 
phase is of first order and there exists sharp phase boundary between spin-liquid and PM 
phases. Anyway, results obtained in this paper support discussions in term of Z 2 models for 
frustrated quantum AF magnets at T = 0. 



IV. SUMMARY 



In this paper we have studied phase structure of the AF spin model on the layered 
anisotropic triangular lattice. We used the Schwinger bosons for represent ting quantum 
spins and also employed the coherent-state path integral methods. We focused on finite- 
T phase diagram and investigated it by means of the MC simulations. We calculated the 
internal energy, specific heat and spin correlation functions. In the absence of the C2-term, 
we found that there exist three phases, i.e., AF, PM and spiral phases. All phase transitions 
between them are of second-order. 

e 
4 

'yi 
2 
: 

0.25 0.5 0.75 1 1.25 1.5 1.75 2 

c 3 




FIG. 20: Encircled region of the phase diagram was studied in the previous paper in terms of 
gauge theories with local Z 2 gauge symmetry [6]. The results obtained in Ref. [6] and the present 
paper show that there exist the PM-tilted-dimer, spiral and deconfined spin-liquid phases in the 
AF magnets on triangular lattice at low T. 

Then we turned on the C2-term and investigated if the deconfined spin-liquid phase ap- 
pears. Calculations of the specific heat and instanton density shows that there exists de- 
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confined spin liquid in the vicinity of the spiral and tilted dimer states. This result is in 
good agreement with the results of our previous study Q in which we assumed a short-range 
spiral order and focused on the region in the phase diagram shown in Figj20l However, the 
present study indicates that there is not sharp phase boundary between the tilted-dimer 
state and spin liquid, i.e., this "transition" is a crossover. At very lower T, it is possible 
that this crossover changes to a genuine phase transition as the imaginary time plays a role 
of another dimension. 

It is very interesting to study how hole doping changes the observed phase diagrams and 
how doped holes behaves in various magnetic phases. This problem is under study and we 
hope that we shall report the results in a future publication. 
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